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Abstract 



The density dependent relativistic hadron field (DDRH) theory is ap- 
phed to strongly asymmetric nuclear matter and finite nuclei far off stabil- 
ity. A new set of in-medium meson-nucleon vertices is derived from Dirac- 
Brueckner Hartree-Fock (DBHF) calculations in asymmetric matter, now ac- 
counting also for the density dependence of isovector coupling constants. The 
scalar-isovector 5 meson is included. Nuclear matter calculations show that 
it is necessary to introduce a momentum correction in the extraction of cou- 
pling constants from the DBHF self-energies in order to reproduce the DBHF 
equation of state by DDRH mean-field calculations. The properties of DDRH 
vertices derived from the Groningen and the Bonn A nucleon-nucleon (NN) 
potentials are compared in nuclear matter calculations and for finite nuclei. 
Relativistic Hartree results for binding energies, charge radii, separation en- 
ergies and shell gaps for the Ni and Sn isotopic chains are presented. Using 
the momentum corrected vertices an overall agreement to data on a level of a 
few percent is obtained. In the accessible range of asymmetries the 6 meson 
contributions to the self-energies are found to be of minor importance but 
asymmetry dependent fiuctuations may occur. 

PACS number(s): 21.65.+f, 21.30.Fe, 21.10.-k, 21.60.-n 
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I. INTRODUCTION 



The modern approach to nuclear structure is based on relativistic models describing 
nuclear matter and finite nuclei as a strongly interacting systems of baryons and mesons. 
Starting from a Lagrangian formulation a phenomenological hadronic field theory is obtained 
by adjusting the meson-nucleon coupling constants to various properties of infinite nuclear 
matter and finite nuclei ll|J^. A connection to free space nucleon-nucleon interactions is 
not attempted. The prototype for such an approach is relativistic mean-field theory (RMF) 
0^ where nuclear forces are obtained from the virtual exchange of mesons, finally leading 
to condensed classically fields produced by nucleonic sources. Using this procedure recent 
RMF models have been remarkably successful in describing nuclei over the entire range of 
the periodic table [0-0] . In order to improve results cubic and quartic self-interactions of the 



mesons fields had to be introduced [pI)| , pA| . Although higher order scalar self-interactions can 
be motivated by vacuum renormalization |^ in practice the strengths of the self-couplings are 
determined in a purely phenomenological way. In mean-field approximation the mesonic self- 
interactions correspond effectively to higher order density dependent contributions. Using 
up to quartic terms a good description of nuclear matter and finite nuclei is obtained but, 
depending on the sign of especially the quartic scalar a self-interaction, instabilities occur 
in the region above saturation density [Q. 

A more fundamental - but also more elaborate - approach is to derive in-medium inter- 
actions microscopically. An appropriate and successful method is Dirac-Brueckner theory 
(DB). Using realistic NN potentials in- medium interactions are derived by a complete re- 
summation of (two-body) ladder diagrams. A break-through was obtained with relativistic 
Brueckner theory which reproduces the empirical saturation properties of nuclear matter 
reasonably well [p!^ -|T8|. Since full-scale DB calculations for finite nuclei are not feasible 
a practical approach is to apply infinite matter DB results in local density approximation 
(LDA) to finite nuclei 0-|2l]]. Retaining a Lagrangian formulation this is achieved by in- 
troducing density dependent meson-nucleon coupling constants taken from DB self-energies 
|T9| . In 10,^ it was pointed out that such an approach does not comply with relativity and 
thermodynamics. A fully covariant and thermodynamically consistent field theory, however, 
is obtained by treating the interaction vertices on the level of the Lagrangian as Lorentz- 
scalar functionals of the field operators. In the density dependent relativistic hadron field 
(DDRH) theory p2| , p3| the medium dependence of the vertices is expressed by functionals of 
the baryon field operators. An important difference to the RMF treatment of non-linearities 
is that the DDRH approach accounts for quantal fluctuations of the baryon flelds even in 
the ground state. Such effects contribute as rearrangement self-energies to the baryon fleld 
equations describing the static polarization of the background medium by a nucleon pB| , P^ . 
Since DDRH theory provides a systematic expansion of interactions in terms of higher or- 
der baryon-baryon correlation functions p2| , pB[ , extensions beyond the mean- level are in 
principle possible. 

In mean-fleld approximation DDRH theory reduces to a Hartree description with density 
dependent coupling constants similar to the initial proposal of Brockmann and Toki |p!9[ . 
The rearrangement contributions signiflcantly improve the binding energies and radii of 
flnite nuclei. Several calculations in the DDRH model for stable nuclei have been performed 
23| , |25| -p7| using density dependent vertices derived from DB calculations with the Bonn A 
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potential p!5|j28[] . Recently, a phenomenological approach to DDRH theory was presented 



^ by determining the density dependence of the vertices empirically. Descriptions of finite 



nuclei of a quality comparable to non-linear RMF models were obtained. Extensions to the 
baryon octet sector and hypernuclei are discussed in []30| . 

The main intention of this paper is to apply DDRH theory to asymmetric matter and 
nuclei far off stability. Since the vertex functionals used in the former applications were 
taken from DB calculations in symmetric matter, information on the density dependence of 
isovector vertices was not available. From the DB results for asymmetric matter, obtained 
in Ref. |]THPT| with the Groningen potential [B^, a new set of coupling constants has been 



derived, now including density dependent isoscalar (a, lu) and isovector (5, p) vertices. 
Contributions from the scalar-isoscalar 6 meson are of special interest at extreme neutron- 
to-proton ratios. 

From inifinite matter calculations it was found that the momentum dependence of self- 
energies has to be taken into account in order to reproduce the underlying DBHF equation 
of state by DDRH calculations. In a strict sense this means to go beyond the static Hartree 
limit. A closer inspection, however, shows that for a mean-field description it is sufficient 
to account for the momentum dependence on an average level. Rather than using the 
self-energies at the Fermi-surface |T^JT^ a more appropriate method is to extract coupling 
constants from self-energies averaged over the Fermi-sphere. This still leads to static vertices 
but incorporating momentum dependent corrections. 

The paper is arranged as follows. In Sec. |T| the DDRH approach and the mean-field 
reduction are reviewed. In Sec. |T| the approach to extract density dependent coupling 
constants from nuclear matter DB self-energies including isovector contributions and the 
momentum corrections is presented. The global properties of the newly determined coupling 
constants are investigated in applications to infinite matter. In Sec. |V| we report on DDRH 
calculations for the isotopic chains of Ni and Sn nuclides, both including stable and exotic 
nuclei. These two isotopic chains are of special interest because several magic neutron 
numbers and the corresponding shell closures are covered on top of the magic Z=28 and Z=50 
proton numbers, respectively. Results for the Groningen and the Bonn A NN potentials 
are compared. Contribution from the scalar-isoscalar 6 meson and the influence of the 
momentum correction are examined in detail. The paper closes in Sec. with a summary 
and conclusions. 



II. DENSITY DEPENDENT HADRON FIELD THEORY FOR ASYMMETRIC 

NUCLEAR MATTER 

A. The model Lagrangian and the equations of motion 

The density dependent relativistic hadron field (DDRH) theory has been presented and 
thoroughly discussed in [^,^,0. In this work we restrict ourselves to a short review of 
the model and to a discussion of the extensions to Ref. [p^ . 

The model Lagrangian includes the baryons represented as Dirac spinors \1/ = {ipp, ipnf ■, 
the isocalar mesons a and cj, the isovector p meson and the photon 7. In addition to former 
models we also include the scalar isovector meson 5 which is important in asymmetric systems 
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and naturally has to be taken into account when extracting the coupling functionals from 
asymmetric nuclear matter DBHF calculations as will be explained in detail later. The 
Lagrangian is 

= Cb + Cm + Cint 
£b = ^[^7m5'^-M]* (1) 

i=a,5 

I E (lF;^:^F^^^'''-ml4^^A^^^') (2) 

e^Q^^'^A^^^^'. (3) 

Here, Cb and Cm are the free baryonic and the free mesonic Lagrangians, respectively, and 
interactions are described by Cint, where 

F^:^ = d,Ai-^ - S^A^f) (4) 

is the field strength tensor of either the vector mesons (k = u, p) or the photon [k = 7) and 
Q is the electric charge operator. 

The main difference to standard QHD models is that the meson-baryon vertices 
Ta (a = a,u!,6,p) are not constant numbers but depend on the baryon field operators \&. 
Relativistic covariance requires that the vertices are functions Ta{p) of Lorenz-scalar bilinear 
forms p(\l/, \1') of the field operators. Two obvious choices are the scalar density dependence 
(SDD) with p = and the vector density dependence (VDD) where fP = j^j^ depends 
on the square of the baryon vector current = \l/7^\l/. In this work we only present results 
for the VDD description since it leads to better results for finite nuclei and gives a more 
natural connection to the parameterization of the DB vertices. This will be discussed in 
detail in the next section. 

As pointed out in |23| the most important difference to RMF or convential DD |19 



theories is the contribution from the rearrangement self-energies to the DDRH baryon field 
equations. This is evident since the variational derivative of C^nt with respect to ^ will also 
act on the vertices. 

~W ~ ~W dp w ^ ^ 

The second term on the right hand side of the equation is the rearrangement contribution 
to the self-energy. Rearrangement accounts physically for static polarization effects in the 



nuclear medium, cancelling certain classes of particle-hole diagrams |24[. The usual self- 
energies are defined as 

S^(°) = f^(p)$. + f5(/5)T$5 (6) 

E'^(o) = f^(p)A(")'^ + Tp{p)TA^P^'' + eQA^^^f". (7) 
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while the vector rearrangement self-energies are obtained from Eq. (^) as 



[ dp 



or, 



dp 
dts 



dp " dp 
Here, is a four velocity with u^u^ = 1. Defining 

the structure of the baryon field equations takes on the standard form 



7m {^d^ 



(m 



^ = 0, 



(8) 



(9) 



(10) 



however, the underlying dynamics is changed by the rearrangement contributions. In addi- 
tion the effective baryon mass M* = M — differs for protons and neutrons due to the 
inclusion of the scalar isovector 5 meson. This is an additional property of the model that 
was not present in the previous formulation [^,^ where only the p meson in the isovector 
part of the interaction was considered. We also find that the density dependence of the p 
and 5 mesons gives an additional contribution to the vector rearrangement self-energies. 



B. Mean-field reduction 

The field equations are solved in the Hartree mean-field approximation. In the Hartree 
approach the highly complex form of the vertex functionals and its derivatives can be treated 



in a simple way using Wick's theorem Calculating the expectation value with respect 



to the ground state |0) the vertices reduce to 

(f,(p)) = TM)) = r.(p) (11) 

and 

/af,(p)\ dT^{p) 



dp / dp 



(12) 



where in the VDD case (p) = p is just the baryon ground state density. 

Meson fields are treated as static classical fields, time reversal symmetry is assumed, 
therefore only the zero component of the vector fields contributes. The meson field equations 
reduce to 

(-V2 + m2)$, = r.(p)/ (13) 

{-V' + mDAt"^ =TMp (14) 

{-V^ + m])^s = Ts{p)pl (15) 

(- + ml)A\^^ = r,(p)p3 (16) 

-V^4"^ = -ep, (17) 
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where the densities are the following ground state expectation values 

P' =_(^^) =Pn+Pl (18) 

p = (^70^) =Pn+Pp (19) 

pl H^rsvl/) =p-:-p; (20) 

P3 = (^70^3^) = Pn-Pp (21) 

and the indices n and p stand for neutrons and protons, respectively. We will use the index 
b = n,p to distinguisch between different nucleons. 

The Dirac equation, separated in isospin, is the only remaining operator field equation 

[7^ {td^ - E^p)) - (M - E^(p))] = (22) 

and contains now the static density dependent self energies 

j:f\p) = r^(p)$. + r,r,(p)<i>, (23) 

^Tip) = r.(p)4^^ + r,r,(p)4^) + ei^4^) (24) 

The self-energies differ for protons and neutrons (r„ = +1, = —1) while the rearrange- 
ment self-energies are independent of the isospin. 



III. RELATIVISTIC HARTREE DESCRIPTION OF INFINITE NUCLEAR 

MATTER 

A. Properties of infinite nuclear matter 

In nuclear matter the field equations further simplify assuming translational invariance 
and neglecting the electromagnetic field. Solutions of the stationary Dirac equation 

[^,k;^-ml]ul{k)=0 (26) 

are the usual plane wave Dirac spinors 

where Xb is a two-component Pauli spinor and the index b distinguishes between neutrons 
and protons. The effective mass ml = M — differs for neutrons and protons due to the 
inclusion of the S meson in the scalar self-energy. The kinetic 4-momenta kl'^ = k^ — H'^ 
and the energy of the particle are related by the in- medium on-shell condition kl^ = m^^ 
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leading to = {kff = ^fef + ml'^. Integrating over all states k < kp^ inside the Fermi 

sphere and introducing Epf, — ^Jkp^ + the scalar and vector densities in infinite nuclear 
matter are found as 



Pb 
Pi 



1.3 

(27r)^ J\k\<kF Stt^ 



-kY J\k 



(28) 



(27r)3 J\k\<k,, El 



27r2 



kp.Ep, + rriu In —2 — 



ml 



(29) 



Calculation of the energy density and the pressure from the energy-momentum tensor 



is straightforward and the results are obtained in closed form 

1 
4 



(30) 



b=n,p 



+ 2 



1 



Y: [^Ep^Pb + mlpl] 

b=n,p 



+ 



b=n,p 



^'^[Pb^'^'+Pt^r 



(31) 



p - ^ h^i = E ] [EF,Pb - mipi] + y: pi>^'^'^ 



i=l 



b=n,p 
1 

2 



b=n,p 



= E 4 [EF.Pb - mlpl] + pE°(^) 

b=n,p 



+ E ^h^r-p^^ 



s(0) 
6 



(32) 



b=n,p 



From these relations it is seen that rearrangement does not afi^ect the energy density but 
contributes explicitly to the pressure p. It is obvious from this that not taking into account 
rearrangement would violate thermodynamical consistency because the mechanical pressure 
p obtained from the energy-momentum tensor must coincide with the thermodynamical 
derivation 



_ 2d (e' 

Pthermo — P T^T \ ~ 
dp \Py 



p. 



(33) 
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B. Nucleon-meson vertices from Dirac-Brueckner theory and the momentum 

correction 



Since the DDRH nucleon-meson vertices are deduced from microscopic Dirac-Brueckner 
calculations the question arises about the best ansatz for the extraction of the DB results. 
This has been discussed extensively in [RO]. For a given infinite nuclear matter DB vertex 



^a{pnm) the mapping to the field theoretical formulation is defined by |23 



Ta{p) = Ta{pnm)S{pl.rn- P^)2pnmdpnm (34) 

Jo 

and directly allows us to apply the DB results to our model. Still, we have not defined 
how to extract the DB vertices Ta{pnm) at a given density pnm from the results obtained 
in Brueckner calculations. Results of Brueckner calculations are the binding energy and 
the DB selfenergies Tp^. The latter are usually calculated by projecting the T matrix 
onto a set of Lorentz invariant amplitudes |[T3|j3^j35|l . They can then be related to coupling 
constants used in mean-field theory as was examined in the local density approximation 



(LDA), e.g. [15^20, 36|. This is usually done on the level of the infinite nuclear matter 
meson field equations by setting the mean- field self- energies equal to S'^^ = T^cpa = 
Plugging this into the meson field equations we find m^S°^ = T^^Pa, with pa being the 
corresponding density to a meson field (j)a as defined in Eqs. ([T8| ) - (|2l|) . 

From Brueckner calculations in asymmetric nuclear matter scalar and vector self- energies 
for protons and neutrons are given, allowing us to extract the intrinsic density dependence 



of isocalar and isovector meson-nucleon vertices. One finds 

(35) 
(36) 
(37) 
(38) 

From this follows that in general the vertices are functions of the Fermi momentum and 
the scalar and vector densities. Specific parameterizations will be discussed in the next 
section. In mean-field theory, only the ratios ^ determine the properties of the EoS. The 
same still holds true for DDRH theory as long as the ratios —{kp) are the same functions 
oi kp Comparing this approach to nonlinear RMF models ||7|, p!0| , p!T| , one finds that 





.2 iK^^^Kkp) 


+ s^(^^)(M 






+ p; 










' 2 Pn 


+ pp 




.2 iK^^^Kkp) 


-^l^'^'^Kkp) 




' 2 


-Pp 




1^ _iEr^)(M 


-^T^\kp) 




/ 2 p„ 


- pp 



the nonlinear o" or u; terms can also be interpreted as density dependent a oi u masses 
or vice versa. However, for finite nuclei this is no longer correct since the rearrangement 
dynamics alters the local single particle properties during the self-consistent calculation. As 
a consequence mass and coupling strength influence the system independently. In this paper 
we only consider constant meson masses and put the medium dependence completely into 
the coupling constants. 
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Self-energies of Brueckner calculations are in general momentum dependent. But the 
usual approach is to neglect the momentum dependence and take the value at the Fermi 
surface |T8| , |31| or to neglect it already a priori |T^. Since the mapping is done on the 
Hartree level, exchange contributions are implicitly parameterized into the direct terms. In 
order to quantify the error from neglecting the momentum dependence, we expand the full 
DB self-energies around the Fermi momentum. 



^DB 



~dk^ 



\k=kf 



+ 0{k' 



(39) 
(40) 



It is common practice to identify the first term with the Hartree self-energy ||T5[. A 
measure of the momentum dependence around the Fermi surface is provided by the second 
term. A quadratic dependence on the momentum has been chosen as supported by Brueckner 
calculations [18|. However, using only TP^{kp) for the determination of the vertices will, in 
general, not reproduce the DB EoS. Up to now a satisfactory solution to this known problem 
20| , |27[| was not yet found. To tackle this problem we introduce as an additional constraint 



that the self- energies have to be chosen such that e^^ = gi^oRH 

As can be seen from Eq. (0), the mean- field contribution of the vector self-energy to the 
potential energy of symmetric nuclear matter is given by ^'^oti^p) = pS°^°^(/cf)- Averaging 
the same contribution from the DB self-energies Tp^{k, kp) over the Fermi sphere and 
requiring it to equal the mean-field potential energy we find the condition 



p(A;^)S°W(fc^) 



(2vr) 



d^kE 



DB 



k\<kF 



{k, kp) 
2 



p{kp)TP''{kp)-T.'{kp)^kl 

15 



p{kp)TP''{k, 



1 



3 



^\kp) 
^TP^{kp 



(41) 



Obviously, the term in brackets is the correction that has to be taken into account in 
order to reproduce the EoS. It should also be included in the extracted self-energies. In 
principle, Ti'{kp) is known from Eq. ( p9D but usually not extracted from DBHF calculations. 
Therefore, we use the approach to calculate Ti'{kp) numerically by adjusting the DDRH 
binding energy to the DB EoS. This can also be interpreted as modifying the vertices 



T\kp) ^ f\kp) ^ T\kp) 



1 



3^^SDB(A;^) 



(42) 



As a first approximation we assume the ratio T,' {kp) /TP^{kp) to depend weakly on kp 
which motivates the introduction of momentum corrected nucleon-meson vertices 



T^{kp)=V^{kF)Jl + Cakl 



(43) 



with Ca being constants determined by adjusting to the DBHF EoS. The rearrangement 
terms are modified as follows 
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+ 1^ r.(M (44) 



A more general ansatz would be to let Co(^f) depend on the Fermi momentum. 

The contribution from the scalar mesons can be treated accordingly, however one has to 
note that the scalar self-energy is also contained in the effective mass m*. Therefore a change 
of E*'-^'' also affects ps and couples back to the modified self-energies. For this reason it is 
not possible to give a closed form for the exact momentum correction. Still, the ansatz from 
equation ( ^31 ) can be used to modify the scalar coupling constants but with the constants 

to be fixed numerically. 

It should be noted that the modified self-energies do not represent the exact DB self- 
energies since the mapping was not done on the single particle level but by adjusting the 
bulk binding energy of the EoS. Also, in order to have momentum dependence on the single 
particle level, exchange terms would have to be taken into account explicitly. However, these 
are already implicitly included in the Hartree vertices through a Fierz transformation P,|5D 



Nevertheless, bulk properties of the DB calculations are retained without adjusting every 
single self-energy separately since C,a was chosen to be a constant thus keeping the number of 
new parameters on a minimal level. The quality of this approximation is measured directly 
by the agreement of the two equations of state. 

One can imagine several possibilities how to calculate the momentum correction. One 
way is to adjust Co- and to the minimum of isospin symmetric infinite nuclear matter. 
Another way would be to keep e.g. the a vertex fixed and adjust Co; for each point of the 
EoS which leads to a density dependent correction Cuiikp). Apparently, the procedure to 
determine the momentum correction is not unique, as already pointed out in P^ , pT| . 



In the next section we are going to present a parameterization of coupling constants 
derived from DB calculations in asymmetric nuclear matter and discuss results of the mo- 
mentum correction by assuming C^a = const. 

C. Fit of the nucleon-nucleon vertices 

Several parameterizations of density dependent coupling constants exist. But they either 
only include density dependence in the isoscalar channel pO[ due to the lack of asymmetric 



nuclear matter DB calculations or they are purely phenomenological P3. Here we present 



a parameterization of asymmetric nuclear matter results |T8|,|3T|] derived from the Groningen 



potential [ll4| , |32| . The mapping of the DB self-energies is done as proposed in equations (pSQ 
- (^) leading to a density dependence in both the isoscalar (a, a;) and the isovector (5, p) 
channel. Figure |l| and Fig. || show the dependence of the ratios —{kp) on the density p(/cf) 
for different asymmetry ratios = pp/ p. In the isoscalar channel the dependence on the 
asymmetry is negligible, in the iscovector channel it is extremely weak, especially around the 
saturation density of about po = 0.16/m~^. We therefore choose the ansatz that the coupling 
constant only depend on the total vector density pi^kp) and not on the proton and neutron 
densities separately. In order to take into account a maximum of information from the DB 
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calculations we fit the average of the self-energies for the asymmetry ratios Og = 0.2, 0.3, 0.4 
Tsj . In |20| a polynomial expansion of in kp around the saturation density kp^^ was 



chosen, leading to an excellent fit of the self-energies around po- But due to the polynomial 
approach the behavior at very low and very high densities is not perfectly stable. Intending 
applications over wide density regions ranging from nuclear halos to neutron star conditions 



in future investigations we choose a rational approximation as proposed in Ref. |29 



r.(p) 



l + c„(^ + e,) 



(45) 



A clear advantage of such a rational form is the well defined behavior at low and high densities 
turning into a constant at very high densities. The results for the fit for po = 0.16/m~^ 
are displayed in Fig. |I| for the isoscalar channel and in Fig. |^ for the isovector channel of 
the interaction. The parameters are shown in Table |. The description of the DB results is 
very good, in the isoscalar channel it is even sufficient to require da = Ca- In the isovector 
channel the fit is more difficult, especially since the (5-meson has an ascending slope at high 
densities. This requires the additional parameter to describe the low and the high density 
behavior equally well. 



D. Results for infinite nuclear matter 

To check the quality of the effective parameterization (Table D of the in-medium de- 
pendence of the vertices it is instructive to look at the infinite nuclear matter EoS. The 
calculation was done by solving the meson field equations ( p!^ - |T6|) with the densities from 
equations (^Sf ) and (p9l). Results for symmetric nuclear matter (a^ = 0.5), pure neutron 
matter (a^ = 0.0) and for nuclear matter with an asymmetry ratio of = 0.3 are shown 
in Fig. ^. Displayed are also the DB binding energies from ||18[ for the same NN potential. 



One sees that the equation of state is clearly not reproduced even though the fit describes 
the self-energies at kp very well. While the DB EoS has a binding energy of e/po = —15.6 
MeV and a saturation density po = 0.182 fm~'^, the standard choice leads to a Hartree EoS 
which is about 2.5 MeV weaker bound (e/po = —13.13 MeV) and the saturation density is 
shifted to lower densities (po = 0.161 fm~^). As discussed in Sec. [ill B| due to the approxi- 



mations made when neglecting the momentum dependence of the self-energies this could be 
expected. 

We apply our momentum correction scheme to the coupling constants in a two step 
process. First we restrict ourselves to symmetric nuclear matter and try to reproduce the DB 
EoS by adjusting and F^. This is done in the constant momentum correction scheme by 
choosing and in such a way that the saturation point of DB calculations is reproduced. 
We find the very small corrections Co- = 0.00804 fm~^ and (uj = 0.00103 fm~^ and are able to 
reproduce the EoS very accurately as can be seen in Fig. ^. It is important to note that even 
though we only adjusted one point we are able to reproduce the binding energies at low as 
well as at high densities. This justifies our assumption of a kp dependence for the correction 
of the coupling constants. Figure ^ compares the momentum corrected couplings to the 
original ones. The correction increases with higher densities (or momenta) as expected from 
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the functional form of the couphng constants but remains small for the a meson and nearly 
negligible for the uo meson. Nevertheless these small corrections suffice to gain 2.5 MeV 
binding energy at saturation density. We conclude that the EoS reacts extremely sensitive 
to small changes in the coupling constants, therefore great care has to be taken when fitting 
the self-energies. In addition the same is appropriate for DB calculations. One has to be very 
careful when extracting the self-energies and needs a consistent scheme for the projection 
onto the Lorentz invariants. 

After fixing To- and the second step is to adjust the couplings in the isovector channel. 
This is done by keeping fixed and adjusting C,p for each given DB binding energy to neutron 
matter. In this approach one obtains a density dependent correction (p{kp). The correction 
is incorporated in the DB self-energies and the p meson-nucleon vertex is readjusted. The 
corrected self-energies and the fit through them are shown in Fig. ^, the parameters are given 
in Table |Il|. From Fig. ^ it is seen that the new fit reproduces the EoS of neutron matter very 
well and this even at high densities where the static fit to the DBHF self-energies leads to an 
interaction being far too repulsive. This is important for applications to neutron stars where 
the high density behavior plays a crucial role. Figure ^ also shows an excellent accordance of 
the calculations with DB results at intermediate asymmetry ratios, e.g. = 0.3, especially 
around the saturation density. This is a very important result since the interaction was only 
adjusted to pure neutron matter and justifies our assumption that the parameterization of 
meson-nucleon vertices is asymmetry independent. We also confirmed this for other values 
of as where DB results were available from [|18| . 



In Table |T| nuclear matter properties for the presented models are given. Saturation 
density and binding energy of the momentum corrected DDRH calculation reproduce the 
DB data very well, also the asymmetry-energy coefficient = 26.1 MeV, determined by 

d e 

is in compliance with DB value of 25 MeV even though this value was not taken into account 
for the adjustment of the isovector interaction. 



IV. RELATIVISTIC HARTREE DESCRIPTION OF FINITE NUCLEI 

A. Properties of finite nuclei 

The density dependent interaction derived in the preceding section for nuclear matter is 
now applied to finite nuclei in Hartree calculations. We solve the full meson field equations 
( p!3| - |T7|) and the Dirac equation (|^) in coordinate space. The Dirac equation is solved for 
the upper and lower components of the eigenspinor xpb simultaneously. The set of coupled 
equations is solved self-consistently under the assumption of spherical symmetry. Pairing 
effects in the particle-particle (p-p) channel have to be taken into account in open shell nuclei. 
Since we are mainly interested in the mean-field particle-hole (p-h) channel and especially 
in the isovector properties of the interaction, the BCS approximation was used. This is a 
standard procedure in relativistic and non-relativistic mean-field approaches. Following P? 



a constant pairing matrix element of G = 2.15 MeV/-\A4 was assumed and the standard 
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set of BCS equations |3§] was solved independently for protons and neutrons, respectively. 
In exotic nuclei close to the dripline pairing effects can be very important if the Fermi 
energy is close to the continuum. Here, relativistic Hartree-Bogoliubov (RHB) calculations, 
taking into account the coupling of bound states to the continuum, have been performed 



for phenomenological interactions [p9|-^ leading to excellent agreement with experimental 
results. We found in our calculations that pairing gives only minor to negligible contributions 
compared to the effects from the microscopic interaction in the p-h channel (less than 2% 
to the total binding energy of the Ni and Sn isotopes). 

The center-of-mass correction which gives a significant contribution to the binding energy 
of light nuclei is treated in the usual harmonic oscillator approximation 

Ecm = ~hu with hu = 41A-^/^MeV. (47) 

Then, the total ground state energy that has to be compared with experimental data is 
given by 

Eg,s. = Emf + -E'pair + Ecm (48) 

where the Hartree ground state energy is obtained from the energy momentum tensor 
through spatial integration of its component. 



i,ei<ep 



b=p,n 

+^Jd'rpp{r)eA^^\r). (49) 

The Si are the Dirac eigenvalues of particles in positive energy eigenstates and energies less 
or equal the Fermi energy ep and the vf are the occupation probabilities obtained from BSC 
pairing. 

In order to examine the effects of the density dependent isovector coupling constants, 
it is instructive to look at the mean-field potentials. In leading non- relativistic order the 
effective central potential is given by the difference of the strongly attractive scalar and 
repulsive vector fields 

f/f = - (50) 
while the strength of the spin-orbit potential 

TjSO ^ 1 —Or (Sfo + 11°) , 

^ 2ME + M- (Sg + SO) ^ ' 

is determined by the sum of these fields. 

Interestingly, the spin-orbit potential differs for protons and neutrons since the self- 
energies depend on the isospin. We define the isoscalar and isovector spin-orbit potentials 

as 
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^o"^ = ^ K + t/r) (52) 

as a measure for isovector spin-orbit interactions. We expect to see an enhancement of the 
isovector potential U^'^ due to the inclusion of the 6 meson. While, to a large extend, the 
contributions of the p and the 5 mesons compensate each other in the central potential U^, 
producing an effective isovector potential that is comparable in strength to the one obtained 
in calculations that include only the p meson, the isovector self-energies of the p and S 
mesons add up and increase the isovector spin-orbit potential U^'~'. This will be discussed 
in detail in the next section. 



B. Closed shell nuclei 



As a first test of the momentum corrected density dependent interaction we examine 
closed shell nuclei. Nuclei considered in our calculations were the doubly magic nuclei '^^O, 
^°Ca, ^^''Sn, ^'^^Sn, ^^'^Pb, having major shell gaps at 8, 20, 50 and 82 for both protons 
and neutrons, as well as nuclei where one or both types of nucleons have only a minor 
(semimagic) shell gap at 28 or 40 (^^Ca, ^^Ni, ^^Ni, ^'^Ni, ^°Zr). The nuclei have measured 



binding energies with very small errors except for the proton rich nuclei. The binding 
energy of ^''^Sn has been measured with a relatively large error '^^Ni, whose binding 
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energy can be extrapolated in terms of the mirror binding energy difference to ^^Ca 
has recently been produced experimentally The set of nuclei includes stable as well as 
(3 unstable nuclei covering the range from the proton dripline to very neutron rich isotopes 
and allows us to investigate the interaction, determined from asymmetric nuclear matter, 
extensively on finite nuclei. 



1. Binding energies and charge radii 

In Table we display results for closed shell nuclei calculated with the DDRH param- 
eterization of the Groningen NN potential derived in Sec. |T|. Figure |^ shows the relative 
error for charge radii and binding energies compared to experimental results. The parame- 
terization derived directly from the self-energies without momentum correction completely 
fails to describe the experimental results. While the charge radii are described very well, 
all nuclei are strongly underbound by about 30%. The proton dripline of the tin isotopes 
is already reached before ^°°Sn due to the weak binding of the protons. But this had to be 
expected from the infinite nuclear matter results where the binding energy at saturation is 
also too weak. Applying the momentum correction from Sec. [Ill D| , the results are greatly 
improved. While the binding energies are still underestimated by about 10%, compared to 
the static parameterization the improvement is remarkable. On the other hand, the size 
of the charge radii decreases, but with an error of about 4% the general agreement with 
experimental data is still satisfying. 

Comparing our results with density dependent calculations for the Bonn A interaction 
PB| , |25|P7| , especially the description of the binding energies is less satisfactory, but the 
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reason for this lies in the different NN potentials. While the parameterization of the Bonn 
A NN potentials has at a nuclear matter saturation density of po = 0.162 fm~^ an energy 
of about —16.3 Mev and for this reason tends to overbind finite nuclei, the saturation point 
of the Groningen NN potential is shifted to higher densities and the saturation energy is 
only about —15.6 MeV. These properties are clearly translated to finite nuclei, leading for 
the Groningen interaction to an underestimation of the binding energy, and a high nucleon 
density and deep mean-field potential inside the nuclei. The high saturation density stronger 
localizes the protons inside the nuclei leading to smaller charge radii. We compared the 
theoretical charge density to experimental results and found an overestimation of about 5% 
to 10% for all considered nuclei. Therefore, to improve results for finite nuclei, Brueckner 
calculations and NN potentials have to be improved first. 

We were also interested in seeing how sensitive results for finite nuclei are to the momen- 
tum correction and if it is possible to improve results by adjusting the momentum correction 
factors (a to finite nuclei. No attempt was made to optimize the parameters in order to ob- 
tain a perfect fit of data. Isovector coupling constants remained unchanged. We found that 
a correction of Co- = 0.008 fm~^ and (^j = —0.002 fm~^ provides a reasonable description 
of both charge radii and binding energies. This modification essentially corresponds to a 
weakening of the repulsion of the vector self-energy compared to the original nuclear matter 
parameterization, but other modifications lead to similar results. It should be noted that 
our discussion of the momentum correction in infinite nuclear matter only partly applies 
to the (a determined for finite nuclei and that the strict relation to the factor S'/S^^ is 
not retained. Obviously, we implicitly take into account higher order effects and especially 
correct for the deficiencies of the NN potential in reproducing the binding energies of finite 
nuclei. This is demonstrated by the reversed sign of compared to the nuclear matter 
value. 

Results are shown in Table 0. As can be seen from Fig. charge radii are further 
decreased while the agreement with the experimental binding energies is satisfying. In 
general, an improvement in the binding energies leads to smaller values for the charge radii 
and it is not possible to reproduce both observables simultaneously with the same Groningen 
parameter set. An improvement of the results would be possible by adjusting the a mass 
to better reproduce the charge radii but a phenomenological fit to experimental results was 
not the goal of our investigations. 

We also investigated if the momentum correction could possibly improve the description 
of finite nuclei for the Bonn nucleon potential. We used the parameterization of the Bonn A 
a and u meson-nucleon vertices from that was also applied to a few closed shell nuclei 
in The p coupling strength is chosen as gp/4:7i = 5.19 with a mass of nip = 770 

MeV. Adjusting the momentum correction factors to the values (a- = —0.0030 fm~^ and 

= —0.0015 fm~^ leads to a good agreement with the experimental binding energies 
(except for some N = Z nuclei where p-n pairing might be important) while preserving the 
good results for the charge radii. Results are presented in Table ^ and Fig. 0. One reahzes 
the same behavior as for the Groningen parameterization. Now, a decrease in the binding 
energies leads to an overestimation of the charge radii. 

The modification of the density dependent vertices is very small. Comparing the original 
fit of Haddad and Weigel |2^ to the momentum corrected vertices we found that, in the 
density range important for finite nuclei, the reproduction of the nuclear matter self-energies 
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is of the same accuracy. This leads to the conclusion that the momentum correction is clearly 
able to improve results and confirms our statement from Sec. Ill that bulk properties are 
very sensitive to small changes in the vertices. 

We also recalculated the equation of state for the Bonn A potential with the modified 
vertices. Comparing to the original parameterization we found the agreement with the 
DB binding energy at saturation density to be slightly improved but not perfect. At low 
densities both parameterizations fail to reproduce the EoS. This seems to be a problem of 
the calculations of Brockmann and Machleidt who assumed a priori momentum-independent 
self-energies and fitted them only to the positive-energy matrix elements of the scattering 
matrix. We tried to adjust the coupling constants with our momentum correction procedure 
but were unable to resolve the inconsistency between the DB self-energies and the DB binding 
energy, e.g. to reproduce the EoS for all densities. This is not surprising considering the 
method of the DB calculations and the fact that the momentum correction is not able to 
change the global properties of the density dependence of the self-energies. 

As can be seen from Figs. ^ and ^ the momentum correction always leads to a shift of 
the relative errors that is nearly identical across all considered nuclei. The reason is that 
the general properties of the interaction were not altered, these are already determined by 
the density dependence of the meson-nucleon vertices. Also the isovector part of the inter- 
action remained unchanged. We conclude that the momentum correction is able to improve 
the extraction of the self-energies from Brueckner calculations preserving the microscopic 
structure of the NN interaction. 



2. Spin- orbit splitting 

As stated before, we expect the effects of including the 6 meson to manifest in the 
isovector spin-orbit potential and in isovector dependent effective masses. In Fig. | effective 
masses for a representative sample of nuclei are displayed. While in the N = Z nucleus ^°Ca 
the difference between proton and neutron masses is induced solely by Coulomb effects and 
is negligible, in the neutron rich nucleus ^^^Sn the proton and neutron effective masses differ 
by about 10% at central density. This enhancement is caused by the isovector contribution 
to the scalar self- energy. The mirror nuclei ^^Ca and ^^Ni nicely demonstrate the effect of 
the scalar isovector density by showing a reversed behavior in the effective masses. We found 
that the momentum correction does not strongly effect the effective masses, only slightly 
decreasing (increasing) them due to a deeper (flatter) core potential. 

The effective masses calculated with the Groningen potential are quite small, leading, 
together with the relatively large self-energies, to a large spin-orbit potential as can be seen 
from Fig. ^ Analogous to the effective masses no strong dependence on the momentum 
correction was found. For comparison we also display the spin-orbit potential of the Bonn 
A and the phenomenological NL3 parameter set The Bonn A spin-orbit potential is 
too small with the peak structure at the nuclear surface shifted to larger radii. We found 
the same behavior for the effective mass, being larger and less localized in the center of 
the nucleus. The properties of the Groningen and NL3 spin-orbit potentials are relatively 
similar, but the spin-orbit splitting of the Groningen potential is too large (about a factor 
two larger than the one obtained with Bonn A) as can be seen from Table [VII[ Compared 
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to the NL3 parameterization that describes the experimental values well, Als is about 25% 
too strong, while for Bonn A it is about 30% too weak. 

Figure |1^ illustrates the isovector dependence of the spin-orbit potential discussed in 
Sec. [IV A| . While U^'^ has a non-negligible strength in neutron-rich nuclei and the mirror 



nuclei ^^Ca and ^^Ni its main contribution is not located at the nuclear surface. One could 
expect to see an effect comparing the spin-orbit splitting of neutrons and protons or mirror 
nuclei but we found no systematic behavior. The isovector spin-orbit potential does not 
seem to manifest itself strongly in the single particle energies and is probably overshadowed 
by the bulk isovector potential. 

C. Ni and Sn isotopes 

The isotopic chains of Ni and Sn are of particular interest for nuclear structure calcu- 
lations because of their proton shell closures at Z=28 (Z=50). They also extend from the 
proton dripline that is found nearby the doubly magic ^^Ni and ^°°Sn nuclei to the already 
j3 unstable neutron- rich doubly magic ^^Ni and ^^^Sn isotopes. This allows us to investigate 
the isovector properties of the derived density dependent interactions and to test the inter- 
actions in regions far off stability. An extensive investigation of these nuclei can be found 
in p9| , ^ , ^ and the references therein. 



1. Binding energies 

Theoretical and experimental binding energies per nucleon are compared in Fig. |TT| for 
nickel and in Fig. |1^ for tin. The Groningen interaction derived from nuclear matter strongly 
underbinds the Ni and Sn isotopes. The reason for this was already discussed in Sec. [IV B| . 



The adjusted parameter set describes neutron rich nuclei reasonably well but fails on the 
neutron-poor side. Agreement for the Sn isotopes is relatively poor. This is mainly caused 
by the extremely strong shell closure at N=82 that shifts the minimum of the binding 
energy from the experimental value of approximately ^^^Sn to the doubly magic ^^^Sn. The 
reason is the too strong spin-orbit splitting of the Groningen potential that affects the 
subshell structure of the Sn isotopes and enhances the shell gap. Nevertheless, the correct 
shell structure is reproduced as can be seen from Fig. where the two-neutron separation 
energies 

S2n{Z, N) = B{Z, N) - B{Z, AT - 2) (54) 

are compared with experimental values. The agreement is satisfying except for the too 
strong shell gaps for ^^^Sn and ^^Ni. The momentum correction does not strongly affect 
S2n since it mainly shifts the binding energy of the isotopic chain not affecting the binding 
energy differences. 

The Bonn A parameter set describes the neutron-poor Ni isotopes fairly well but 
overbinds the neutron-rich Ni and the Sn isotopes. On the other hand, agreement of the 
results for the momentum corrected interaction with experimental data on the neutron- rich 
side is very good but binding energies of neutron-poor isotopes are underestimated by 0.2 



17 



MeV. Also, comparing the two-neutron separation energies (Fig. p^ ), experimental results 
are reproduced well, but some deviations can be noticed. The shell gaps for ^^^Sn and ^^Ni 
are relatively small whereas the gap for ^^Ni is too large. This also explains why the min- 
imum of the binding energy is found at ^^Ni and the isotopes of the ^^Ni to ^^Ni shell are 
underbound. The reason for this probably lies in the weak spin-orbit splitting of the Bonn 
A potential. 

Realizing that for N ^ Z nuclei only the Coulomb interaction determines the difference 
between the neutron and proton energies while the isovector interaction is strongly sup- 
pressed we conclude that the underbinding of the neutron- poor nuclei is caused by a too 
weak central isoscalar potential. On the other hand the neutron- rich isotopes are mostly 
affected by the isovector interaction that acts attractive for protons, balancing the Coulomb 
repulsion, and repulsive for neutrons. Since neutron-rich isotopes are described well this 
is an indication that for both the Groningen and the Bonn A potential the microscopic 
isovector potentials seem to be too weak. One can visualize this easily by shifting the curves 
in Figs. |ll] or |l^ vertically to reproduce the binding energy of N = Z nuclei. This in- 
creases in first order only the isoscalar central potential. Then, the neutron-rich isotopes 
are overbound, indicating that the repulsion from the isovector potential is too weak. 



2. Density distributions 

Self-consistent neutron density distributions and charge densities for the Ni and Sn iso- 
topes are displayed in Fig. |14[ Closed shell nuclei are marked by bold lines. The distributions 
are only presented for the Bonn A potential but for the Groningen potential the same sys- 
tematics is obtained. For the calculation of the charge density the theoretical point particle 
density distribution pp is folded with a Gaussian proton form factor |^ with yij^p = 0.8 
fm. From ^'^Ni to ''^Ni the charge density in the interior is reduced by about 30% accompanied 
by a mild increase of the charge radius by about 5%. For the neutron densities a more drastic 
evolution is found. Beyond ''^Ni (N=50 shell closure) the 2^5/2 subshell is being filled and a 
thick neutron skin is build up. This leads to sudden jump in the neutron rms radii while the 



proton rms radii increase only slowly. This behavior is more clearly visible in Fig. |T5| where 
the differences of the proton and neutron rms radii are shown. Approaching the proton 
dripline a relative thick proton skin below N=28 is predicted. We also compare our results 



to calculations with the NL3 interaction. In Ref. |^ relativistic Hartree-Bogoliubov calcu- 



lations for the Ni and Sn isotopes using this interaction were performed leading to excellent 



agreement with experimental data. Figure 15 shows very good agreement for the Bonn A 



potential with the results for the NL3 interaction. The Groningen potential has the same 
tendency but with smaller values. This is explained by the systematic underestimation of 
the rms radii (Fig. |^) that also causes a reduced difference between proton and neutron rms 
radii. Nevertheless, the same neutron skin is found. Since the Groningen potential leads to 
larger shell gaps for reasons of its strong spin-orbit potential, at the (sub)shell closure N=40 



a small increase is already seen. Figure |T5| also illustrates that the results are independent 
of the momentum correction. This is obvious since the correction only shifts all rms radii 
to higher or lower values. 
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The same calculations were performed for the Sn isotopes. In Fig. |14| density distribu- 
tions from ^'^"Sn to ^"^"Sn are displayed showing the same neutron skin as the Ni isotopes. 
Here, one finds a sudden jump beyond N=82 where the l/iii/2 shell is filled and the 2/7/2 
subshell becomes populated (Fig. ^). Again the agreement between the Bonn A and the 



NL3 parameter set is excellent. In Ref. calculations for the Sn isotopes with non- 



relativistic interactions were performed and identical results were found. This shows that 
the observed neutron skins in Sn and Ni are relatively independent of the NN interaction 
and the theoretical approach. 



V. SUMMARY AND CONCLUSION 

We have extracted density dependent meson-nucleon vertices from DBHF self- energies 
of asymmetric nuclear matter, derived from realistic NN potentials. For the Groningen NN 
potential we found that the coupling constants can be expressed solely in terms of the nuclear 
matter baryon density p and are independent of the asymmetry fraction for which the self- 
energies were calculated. The extraction of momentum- independent self-energies introduces 
errors in the density dependence of the interaction and leads to difficulties to reproduce the 
DB EoS in the DDRH approach. We introduced a momentum correction to account for 
this error and found excellent agreement of our calculated EoS for all asymmetry fractions 
Qs = Pp/p with the Brueckner results. While the corrections are very small, the sensitivity 
of the saturation point of nuclear matter and the binding energies of finite nuclei to this 
correction is very high. This indicates the difficulty of extracting momentum independent 
coupling constants. 

Applying the momentum corrected interaction to finite nuclei, we find that the agreement 
with experimental results is satisfactory, taking into account that the used parameterization 
contains no free parameters for finite nuclei. For the Groningen parameter set, the binding 
energy per nucleon is underestimated by 0.5 to 1 MeV whereas the charge radii are about 0.1 
to 0.15 fm too small. These results are comparable with other mean-field calculations p5|-27[ . 



Using the Bonn A potential results are improved, indicating that the Groningen potential is 
not attractive enough, especially in the low density regime. We also adjusted the momentum 
correction factors to improve the description of the properties of finite nuclei. This leads to 
a good agreement with experimental data while keeping the excellent reproduction of the 
DB self- energies. 

The inclusion of the 5 meson introduces different effective masses for protons and neu- 
trons and strongly enhances the isovector spin-orbit potential. However, we found no sys- 
tematic effect in the isovector spin-orbit splitting. In order to further examine this effect, 
detailed experimental data for the single- particle energy levels and for exotic nuclei are 
necessary. Compared to phenomenological RMF interactions, we found the spin-orbit split- 
ting of the Groningen potential to be enhanced and the one of the Bonn A potential to be 
suppressed. This becomes visible in the shell structure of some exotic nuclei. We also found 
indications that the isovector interaction of both microscopic interactions seems to be too 
weak to reproduce the complete isotopic chains of Sn and Ni correctly. 

In general, we think that the results of the DDRH approach are quite satisfactory and 
that the momentum correction provides a consistent scheme to reproduce DBHF calculations 
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and improve the agreement with finite nuclei. Improvements of the results could possibly 
be achieved by going beyond the ladder approximation and including, e.g., three-body in- 
teractions and ring diagrams. In future investigations we also plan to apply the density 
dependent interactions to neutron stars to gain additional insights in the properties of the 
isovector density dependence. 
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200 I 1 1 1 1 1 1 1 1 1 1 1 1 i > r-^ I 

0.0 0.1 0.2 0.3 0.4 

P [fm '] 

FIG. 1. Effective density dependence of the a (upper part) and the lo (lower part) me- 
son- nucleon vertices. Shown are results extracted from DB selfenergies from the Groningen NN 
potential ||l^ calculated for the asymmetry ratios = 0.2, 0.3, 0.4. The solid line is the asymmetry 
independent fit through the average of the DB results (open squares). 
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FIG. 2. Same as Fig. [l| but for the p (upper part) and 5 (lower part) meson-nucleon vertices. 




-20 I I 
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nuclear density p [fm"^] 

FIG. 3. Equation of state for different asymmetry ratios calculated from the Groningen NN 
potential. DB results from are represented by symbols. The dashed line shows the DDRH 
EoS for meson-nucleon vertices fitted to the DB self- energies, the solid line denotes results with 
momentum correction. 
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FIG. 4. Momentum corrected density dependence of the a and u meson-nucleon vertices. Re- 
sults are shown for = 0.00804 fm~^ and Co.) = 0.00103 fm~^ (solid line) and compared to the fits 
from Fig. |l] (dashed line) . 
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FIG. 5. Momentum corrected density dependence of the p meson- nucleon vertex. Results are 
shown for the fit (solid line) through the corrected self-energies and compared to the fit from Fig. ^ 
(dashed line). 
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FIG. 6. Relative errors for charge radii {pc — pf^^)/p'^^ (upper part) and binding energies 
(Eb — E'^'^) / E'^'^ (lower part) obtained with the Groningen parameterization. Shown are results 
for magic and semimagic nuclei. Results without momentum correction are denoted by upper 
triangles, results for ("o- = 0.00804 fm^^ and Co; = 0.00103 fm~^ (adjusted to nuclear matter) by 
squares and results for Co- = 0.008 fm^^ and Coj = —0.002 fm~^ (adjusted to finitite nuclei) by 
circles. The lines are drawn to guide the eye. 
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FIG. 7. Same as Fig. |6| but for the Bonn A parameterization from j2C,23|. Results without 
momentum correction are denoted by squares and results for Ca = -0.003 fm^^ and = -0.0015 
fm~-^ by circles. The lines are drawn to guide the eye. 
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FIG. 8. Isospin dependence of the effective masses of closed-shell nuclei from relativistic density 
dependent Hartree calculations using the parameterization of the Groningen NN potential. 
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FIG. 9. Spin-orbit potentials of closed-shell nuclei from DDRH calculations. Shown are the 
isoscalar potentials Uq'^ for the Groningen parameterization (solid lines) and the Bonn A param- 
eterization (dotted line). For comparison results calculated with the phenomenological NL3 RMF 
parameter set are shown (dashed lines). 
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FIG. 11. Binding energies per nucleon for the nickel isotopes. Shown are results with (open 
circles) and without momentum correction (solid circles) for the parameterizations of the Bonn 
A potential. Results derived from the Groningen parameter set are denoted by solid triangles 
(momentum corrections adjusted to nuclear matter EoS) and open triangles (adjusted to finite 
nuclei). Experimental binding energies are denoted by crosses (taken from Ref. |^^). The lines are 
drawn to guide the eye. 
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FIG. 12. Same as Fig. 11 but showing binding energies of the tin isotopes. 
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FIG. 13. Two neutron separation energies for the nickel (lower part) and the tin (upper part) 
isotopes. Results for the Groningen (triangles) and Bonn A (circles) parameterization are shown. 
Experimental separation energies are denoted by crosses (taken from Ref. [Q). The lines are 
drawn to guide the eye. 
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FIG. 14. Charge and neutron density distributions for the isotopic chains of the Sn and Ni 
nuclei, calculated with the Bonn A parameter set. 
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FIG. 15. Difference of the root mean square radii of neutrons and protons of the nickel (lower 
part) and the tin (upper part) isotopes calculated with the Groningen (triangles) and Bonn A 
(circles) parameterizations. For comparison results calculated with the phenomenological NL3 
RMF parameter set are shown (crosses). The lines are drawn to guide the eye. 
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TABLES 





C 


UJ 


X 



P 


ma[MeY] 


550 


783 


983 


770 


da 


13.1334 


15.1640 


19.1023 


12.8373 


ba 


0.4258 


0.3474 


1.3653 


2.4822 


Cq, 


0.6578 


0.5152 


2.3054 


5.8681 


da 


0.7914 


0.5989 


0.0693 


0.3671 




0.7914 


0.5989 


0.5388 


0.3598 



po = 0.16 [fm~ 



TABLE L Parameterization of the density dependent couplings from equation (|45| ) extracted 



from DB calculations in asymmetric nuclear matter |18 



Oq, 6q, Cqi da da 

19.6270 1.7566 8.5541 0.7783 0.5746 

PO = 0.16 

TABLE IL Parameterization of the p meson-nucleon vertex after adjusting the neutron matter 
DDRH equation of state to the DB binding energy. 



DBHF DDRH DDRH corr. 
Psat [fm-^] 0.182 0.161 0.180 

e/p [MeV] -15.5 -13.13 -15.60 

K [MeV] - 211 282 

m*/M - 0.592 0.554 

04 [MeV] 25 28.2 26.1 

TABLE HL Comparison of infinite nuclear matter properties obtained in the DDRH model with 
results from DB calculations for the Groningen NN potential. Results are shown for calculations 
with and without momentum correction. 
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exp. 




i?/74[MeV/AJ 




2.74 


7.98 


Aft ^ 


3.48 


8.55 


^^Ca 


3.47 


8.67 




4.27 


8.71 


208pb 


5.50 


7.87 


48Ni 




7.27 


56^1 




8.64 


68Ni 




8.68 


lOOSn 




8.26 


i32Sn 




8.26 



TABLE IV. Experimental values of the rms charge radii and binding energies per nucleon for 



(semi)magic nuclei. Data are taken from Refs. [37,42,43,46] 



roningen 




= 0.0 




0.00804 fm-2 




0.008 fm-2 






= 0.0 


C = 


0.00103 fm-2 


Cu = 


-0.002 fm-2 






E/A 




E/A 




E/A 


160 


2.76 


5.65 


2.61 


7.00 


2.56 


7.56 


40Ca 


3.47 


5.96 


3.30 


7.52 


3.24 


8.16 


48Ca 


3.48 


6.09 


3.31 


7.79 


3.25 


8.49 


90Zr 


4.26 


6.00 


4.06 


7.77 


3.99 


8.49 


208p]^ 


5.49 


5.28 


5.24 


7.10 


5.15 


7.84 


48Ni 


3.77 


4.68 


3.54 


6.28 


3.46 


6.95 


s^Ni 


3.73 


5.78 


3.53 


7.60 


3.46 


8.35 


68Ni 


3.88 


6.13 


3.69 


7.83 


3.63 


8.52 


loosn 


(4.46) 


(5.42) 


4.24 


7.25 


4.16 


8.00 


i32Sn 


4.69 


5.86 


4.47 


7.68 


4.40 


8.42 



TABLE V. Root-mean-square charge radii rc[fm] and binding energies per nucleon Eb/A 
[MeV/A] of (semi)magic nuclei from density dependent relativistic Hartree (DDRH) calculations 
using the Groningen NN potential. Results for different momentum correction factors (a and Ct^ 
(see text) are shown. The values in parenthesis correspond to unbound nuclei. 
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Bonn A 


C<7 


= 0.0 




Ca = - 


-0.0030 im 


-2 






= 0.0 




Clo = - 


-0.0015 fm 


-2 








E/A 






E/A 




2.75 




8.58 


2.79 




8.30 


40 /-I 


3.46 




9.02 


3.50 




8.69 


48 /~i 


3.49 




8.96 


3.53 




8.63 


yu2j. 


4.26 




8.98 


4.31 




8.63 


208p|jj 


5.53 




8.17 


5.60 




7.82 


48Ni 


3.84 




7.62 


3.89 




7.31 




3.79 




8.66 


3.84 




8.31 


1 M i 


3.88 




9.05 


3.93 




8 71 


loogn 


4.51 




8.40 


4.57 




8.04 


Sn 


4.74 




8.67 


4.78 




8.32 


TART.F, VT 


Same as Table |V] but with results for the Bonn A NN potential. 






160 




40 Ca 


48 Ca 




48Ni 


Gron. 


8.0/7.8 




8.3/8.0 


8.0/7.7 




7.5/7.4 


Gron. adj. 


9.0/8.7 




9.6/8.8 


8.8/8.5 




8.3/8.2 


Bonn A 


4.2/4.2 




4.6/4.6 


4.0/4.1 




3.9/3.8 


Bonn A adj. 


4.0/4.0 




4.4/4.3 


3.7/3.8 




3.6/3.6 


exp. 


6.1/6.3 




6.3/7.2 


5.6/4.3 







TABLE VII. Neutron and proton spin-orbit splitting Ai5(n,p) in MeV for the Ip shell (^^O) 
and the Id (^'^Ca, ^^Ca and '^^Wi) shell. Results are shown for the Groningen and the Bonn A NN 
potential, for the interaction derived from nuclear matter and for the adjustment to finite nuclei, 
respectively. For experimental values see Ref. 
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